Networks with desired correlations 
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We discuss a simple method of constructing correlated random networks, which was recently 
proposed by M. Bogufia and R. Pastor-Satorras (cond-mat/0306072). The result of this construction 
procedure is a sparse network whose degree-degree distribution asymptotically approaches a given 
function at large degrees. We argue that this convergence is possible if the desired function is 
sufficiently slowly decreasing. 
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Recently, Marian Bogufia and Romualdo Pastor- 
Satorras proposed a natural way of constructing corre- 
lated networks with desired correlations of the degrees 
of the nearest neighbors [1]. This practical algorithm is 
based on the idea that previously was used in Refs. [2-4] 
for constructing uncorrelated graphs and in Ref. [5] for 
building correlated networks (see also Ref. [6]). In all 
these constructions, (i) some weights (fitnesses, desired 
degrees, etc.) are ascribed to vertices, and (ii) each pair 
of vertices is connected with probability which depend 
on these weights. 

In the first version of this work, I reproduced this al- 
gorithm without knowing that it had already been pro- 
posed. So, in the present, corrected and shortened and 
more methodical version I mostly discuss the range of 
validity of the algorithm and features of networks, gen- 
erated by this method. 

The correlations between degrees of the nearest neigh- 
bors in a graph are naturally described by the joint dis- 
tribution of the degrees of end vertices of an edge of the 
graph, P(fc,fc'), P{k,k') = P{k',k), Y.k,k' Pi^.k') = 1. 
The joint degree-degree distribution determines the de- 
gree distribution P{k) of a network: 
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and (fc") =A:Efc,fc' k''-^P{k,k'). 

The algorithm of Boguha and Pastor-Satorras gener- 
ates sparse random networks with desired degree-degree 
correlations. Suppose one wishes to obtain an ensemble 
of graphs with a desired joint distribution of the degrees 
of the nearest neighbors, P{q,q'). One should assume 
that P{q, q') decreases with q and q' sufficiently slowly. 
Let the number N of vertices in each graph of the en- 
semble be large and fixed. 

The procedure [1] is as follows: 



(i) Create N vertices with a sequence of weights {qi}, 
i = 1, . . . , N independently sampled from the distribu- 
tion 



(3) 



where q = [T,q,q' Pi<l,<l')/q 



(ii) Put a link between i and j vertices with probability 
[7] 



Pi] ^piqi,qj) = 



q P{qi,qj) 
NP{q,)P{qj) 



(4) 



In this model, (i) the average degree of an i-th ver- 
tex coincides with q^, and (ii) the distribution of degrees 
of an individual vertex is a relatively narrow function at 
large degrees. This can be proved by calculating the de- 
gree distribution of a vertex with weight q (see Ref. [1]). 
Here, alternatively, we simply find the first two moments, 
ki and {kf) for the degree of an i-th vertex. 

The statistical weights of graphs g G G in the result- 
ing ensemble, written in terms of pij (adjacency matrix 
elements are a^- = 0, 1), are of a rather standard form 



n(g) = n({p,,},{a,,})(xn 
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(compare with the statistical weights of classical random 
graphs and Ref. [6]). Here, = Py/(1 — Pij)- The 
partition function of the ensemble is 



z{G) = z{{p,,})= J2 n(g) = Enn(fe. 



oc 



geG 



t,] 



-' = 1[(^+Sl]). 



(6) 



So, 
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etc. [The averages are over the statistical ensemble: 
{X{g)) = ZiG)-'j:g^^n{g)X{g).] Consequently, ah 
the moments are equal: (a^) = Pij, n > 1. This, in fact, 
is clear, since a^j takes only two values and 1 with the 
probabilities 1 — pij and pij , respectively. Then, 

fci = ( X] ) = X] Pij = ^ X P(l)P(9i ,Q) = Qi (8) 

^ j ' j Q 

[Eq. (4) was used], so the mean degree in the net k = q, 
and 

(kf) = ( (e «^.)' ) = Epij + (Ep ' - E4 • 

(9) 

Note that Eqs. (8) and (9) may be obtained by using 
more naive arguments. Suppose for brevity that a net is 
of three vertices, 0, 1, and 2. Then fco = 1 • [poi(l— P02) + 
(l-Poi)Po2]+2-poiPo2 =Poi+Po2 [compare with Eq. (8)] 
and (fc^) = l2-[poi(l-Po2) + (l-poi)Po2]+22-poiPo2 = 
Poi +P02 + 2poiPo2 [compare with Eq. (9)]. 

Using the expression (4), one can see that in a sparse 
network, the last term on the right-hand side of Eq. (9) 
is of the order of N^^ and so is negligible. Consequently, 
in the thermodynamic limit (A'^ — » 00), 
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(10) 



That is, the distribution of degrees of an individual ver- 
tex i is indeed relatively narrow, if ki is large, and peaked 
at ql. One should take into account the fact that the 
number of edges in the net, which connect vertices with 
weights q and q' , is 



We ask two (related) questions, (i) How should a 
function f{x) behave at large x to satisfy the condi- 
tion: f{x + c^/x)/ f{x) 1 as X — > 00? (ii) How 
should a function f{x) behave at large x to guarantee 
that [f{x + cy/x) + f{x - cVi)]/[2/(x)] ^ 1 as x ^ 00? 
Here, c = const. 

At first sight, the second, more symmetric condition 
may demand less strong restriction on the behavior of 
f{x). However, both the questions have the same answer. 
Indeed, the second and the third, and the higher terms in 
the series f{x + C\/x)/f{x) = 1 -|- c^/x[df{x)/dx]/f{x) + 



^x[<ff(x)/dx^]l[2f{x)] - 



approach zero at large x if 



the same condition is satisfied: f{x) must decrease with 
x slower than e~^. So, if the given function P(k,k') 

decreases slower than, say, e~^~^^, then Prosuit(^, fc') 
asymptotically approaches P{k, k') at large k and k' [8]. 

Moreover, this condition also guarantees that the re- 
sulting degree distribution Prcsuit(fc) asymptotically ap- 
proaches P{k) at large degrees. On the other hand, in 
the region of small degrees, PresM\t{k,k') deviates from 
desired P{k.k'). and Prosuit(^) deviates from P{k). 

Similar arguments are also valid for other network 
constructions of this type [2-6]: desired distributions 
must decrease sufficiently slowly. On the one hand, this 
implies a serious restriction on the range of degree (or 
degree- degree) distributions which can be reproduced in 
such a way. On the other hand, it is the slowly decreas- 
ing distributions that are most interesting. 
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L{q,q') = -N^P{q)P{q'){a{q,q')) = 
^NqP{q,q') = LP{q,q'), 



(11) 



where L is the total number of edges [note the symme- 
try factor 1/2 in the second term of Eq. (11)]. Then 
one can finally conclude that the resulting degree-degree 
distribution of the net is expressed in terms of a given 
function: 



(12) 



Presnltik, k') = P{k + 6{k), k' + S' {k')) , 



where the deviations, \6{k)\ < Vk, \S'{k')\ < Vk', 
are relatively small at large degrees. The relation (12) 
shows that at large degrees the algorithm of Bogufia and 
Pastor-Satorras provides the degree-degree distribution 
Presuit{k,k') = P{k,k'), only if P{k,k') decreases suf- 
ficiently slowly. The following arguments explain what 
this does mean. 
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